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Abstract 

In theories with low energy supersymmetry breaking, the effective potential for squarks 
and sleptons has generically nearly flat directions, V(<p) ~ M 4 (log(0/M)) n . This guar- 
antees the existence of stable non-topological solitons, Q-balls, that carry large baryon 
number, B ^> (M/m p ) 4 , where m p is the proton mass. We study the behaviour of 
these objects in a high temperature plasma. We show that in an infinitely extended 
system with a finite density of the baryon charge, the equilibrium state is not homo- 
geneous and contains Q-balls at any temperature. In a system with a finite volume, 
Q-balls evaporate at a volume dependent temperature. In the cosmological context, we 
formulate the conditions under which Q-balls, produced in the Early Universe, survive 
till the present time. Finally, we estimate the baryon to cold dark matter ratio in a 
cosmological scenario in which Q-balls are responsible for both the net baryon num- 
ber of the Universe and its dark matter. We find out naturally the correct orders of 
magnitude for M ~ 1 ... 10 TeV: r] ~ 10- 10 (M/TeV)- 2 ( J B/10 26 )- 1 / 2 . 
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1 Introduction 



Take a theory containing scalar fields which carry some unbroken global U(l) charge Q. 
Assume that all charged particles are massive. Suppose that the effective potential 
for the U(l) charged scalar fields has a flat direction, V(<f>) ~ M 4 at large 0, up to 
possible logarithmic terms. Then this theory contains absolutely stable non-topological 
solitons, Q-balls, that have a non-zero value of the global charge Q [|IJ. Flatness of the 
effective potential at large is essential for this statement. For Q > 1, the mass of 
the soliton grows as M Q ~ MQ 3/4 and becomes smaller than the mass of a collection 
of free separated particles with the same charge, mQ 1 , independent of the relationship 
between m and M. For theories without flat directions, where V(<f>) grows at large 
as n with n > 2, Q-ball solutions constructed from scalar fields can exist as well [0, |3|, 
but their energy scales with Q as AQ 1 , so that the question of their stability depends 
on the mass of the lightest particle that carries the charge and on the (computable) 
coefficient A. 

A phenomenologically interesting example of a system with these properties is pro- 
vided by supersymmetric extensions of the Standard Model. Here the role of the U(l) 
charge is played by the baryon number, and the role of scalar fields by squarks. The 
MSSM with unbroken supersymmetry has a lot of gauge invariant flat directions with 
V(4>) = carrying baryon number Q. If supersymmetry is broken at a low energy 
scale M, as in gauge-mediated scenarios || , then the potential is lifted at > M in a 
way that V((f>) stays constant || [7]]. Thus, this theory contains stable non-topological 
solitons which carry baryon number ||. The generic MSSM with other types of SUSY 
breaking also contains a lot of Q-ball solutions || [L0|, but with Mq ~ MQ 1 . 

Stable Q-balls in supersymmetric theories have a number of interesting properties. 
Since they represent the most economic way of packing the baryon number, they like to 
absorb the ordinary matter in the form of protons and neutrons and to convert baryon 
number from fermionic baryons to bosonic baryons (squarks) ||11|| . Large stable Q-balls 
can be naturally produced || in the Early Universe via the decay of an Affleck-Dine 
(AD) condensate [|12|] and can contribute to the cold dark matter ||. Experimental 



signatures of relic dark matter Q-balls are spectacular [13 



At zero temperature, Q-balls are the states with minimum energy at a fixed value 
of the charge Q. They do not decay if the mass to charge ratio Mq/Q for them is 
smaller than that for a free particle. At non-zero temperatures all energy levels of the 
system are populated with non-zero probabilities and, thus, the charge can leak out 
from a Q-ball to the outer space, or, equally, can be absorbed from the surrounding 
plasma. Absorption of the charge into a Q-ball decreases the energy of the system, 
while releasing the charge from it increases the entropy. So, the issue of Q-ball stability 
gets more complicated than at zero temperature. 

The aim of this paper is to consider Q-ball thermodynamics in theories with flat 
directions in the scalar potential. We will see that in these models systems with a 



1 



finite density of the baryon number % do not have a standard thermodynamical limit. 
For large enough volumes the initially homogeneous distribution of % 7^ is unstable 
against Q-ball formation at any temperature. On the contrary, if the volume of the 
system is finite and its charge is fixed, Q-balls cease to exist (evaporate completely) at 
temperatures T > Tq, where Tq is volume dependent. 

In the cosmological context, a mere thermodynamical consideration is not enough 
as the system has a finite time scale given by the expansion rate of the Universe, and 
the (non-equilibrium) rate of Q-ball evaporation is essential. We make an estimate of 
this rate and determine under which conditions Q-balls created in the Early Universe 
survive till present. If the only source of baryon asymmetry of the Universe is the AD 
mechanism, then the baryon asymmetry and cold dark matter may share the same 
origin ||. In this scenario, both matter and dark matter are baryonic, though dark 
matter is constructed from squarks packed inside Q-balls. We show that in this case, 
the baryon to cold dark matter ratio appears to be related to very few parameters such 
as the charge of an individual Q-ball, the Planck scale, the SUSY breaking scale M 
and the proton mass. 

Initially, Q-balls were found and studied in different models in Refs. The 
generation, evolution, evaporation and implications of Q-balls at the high tempera- 
tures of the Early Universe were discussed in various models in [|L4|-||25|. In particular, 
in |T9| , |20| it was found that at high temperatures (but still T <C M) Q-balls tend 



to evaporate whereas at low temperatures, they are stable and tend to grow. The 
influence of SUSY Q-balls on phase transitions was discussed in [^3] . In theories where 
supersymmetry breaking comes from the supergravity hidden sector, Q-balls are unsta- 
ble at zero temperature. Nevertheless, they can be produced in the Early Universe and 



are important for baryogenesis [24]. Moreover, they may provide a natural explanation 



for the ratio between neutralino dark matter and the baryonic matter in the Universe 



25] . However, theories with flat potentials, specific for theories with low energy su- 



persymmetry breaking have, to our knowledge, never been analyzed in detail for all 
temperatures up to T > M. 

This paper is organized as follows. Essentially, it follows two distinct lines: a qualita- 
tive discussion of the physically interesting case of the MSSM with order of magnitude 
estimates, and a more formal and quantitative analysis, but mostly applied to a specific 
toy model. In Sec. |2|, we formulate in general terms the main problems considered. In 
Sec. |3|, we discuss the solutions to these problems on a qualitative level. In Sec. f| we 
perform a more quantitative analysis, but in a simple renormalizable supersymmetric 
toy model with flat directions in the effective potential. A number of claims made in 
Sec. H are substantiated in Sec. |[ Finally, in Sec. || we consider the applications of 
the results to cosmology in the realistic case, again on an order of magnitude level. In 
particular, we consider here the stability of Q-balls at high temperatures, their evapo- 
ration rate, and their survival in the expanding Universe. We also estimate the baryon 
to cold dark matter ratio. Some details related to Sec. |] are in the two Appendices. 



2 



2 The formulation of the problem 



Let us consider a system with some conserved global charge Q, such as B or L (baryon 
number or lepton number; the B + L anomaly is not important for this discussion). The 
symmetry which corresponds to this charge is assumed to be unbroken. On a formal 
level, the main questions we are going to address are as follows: 

1. Assuming that the thermodynamical ground state of the system is homogeneous, 
is the conserved charge carried by individual particles in a plasma, or by a scalar 
condensate (Q- matter), in analogy with Bose-Einstein condensation? 

2. Is such a homogeneous thermodynamical ground state stable against small in- 
homogeneous perturbations? If it is, is it the true ground state or a only metastable 
minimum? 

3. In case a homogeneous ground state is unstable or metastable, what is the global 
minimum of the free energy? In particular, under which conditions can Q-balls consti- 
tute such a ground state? 

4. If one starts from an unstable or metastable solution, at which rate does one 
approach the equilibrium ground state? In particular, at which rate do Q-balls form if 
they are stable, or evaporate if they are unstable? 

In order to answer these questions, let us first fix the basic definitions used. The 
thermodynamics of a system with a conserved charge (the case of zero temperature is 
incorporated as a special case) is determined by the grand canonical partition function 

Zfl = e -^ Ty ^ = e? p ™ = Tr e -^-^\ (2.1) 

where 



g / -9.= q= m&. (2.2) 



V V <9/i 

If Q is kept fixed instead of /x, then the thermodynamical potential to be minimized is 
obtained with a Legendre transformation: 

F(T, V, Q) = Vf(T, q) = Q + ^Q = V[-p + m ]. (2.3) 



The relation inverse to eq. Q2.2 ) is then 



" - {2A) 



Given the Lagrangian of the system, one can by standard methods |26j derive a 
Euclidian path integral expression for Z^. In a field theory containing the fields $, it 
is convenient to express the result in terms of an effective action in the usual way. It 
then follows that 

p(T,pt) = -((3V)- 1 r^lln. w/w=0 . (2.5) 
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The relevant extremum is the minimum, as follows from the general thermodynamical 
principle that, for fixed T, /i, the other variables take such values that the corresponding 
free energy is minimized. 

In the case of a homogeneous extremum, the problem reduces to an effective poten- 
tial. Then 

P(T^) = -V,m dv ^ =0 . (2.6) 
The Legendre transform is V q ($) = V^($) + fiq where q = — <9 M V^($), and 

dV q /d<S>=0 ■ 

(2.7) 

With these tools, the first three of the questions formulated above can, in principle, be 
answered. The first of the questions is particularly straightforward: one just computes 
the effective potential for a fixed charge, V q (&), and inspects whether a non-zero value 
for $ serves to minimize it or not. If yes, then one can have a Q-matter condensate. 

To address the second question, one can study the standard thermodynamical sta- 
bility conditions: these follow from the requirement that the entropy of any subsystem 
must be at a maximum with respect to local fluctuations of temperature, volume, and 
the number of charged particles. Applied to the present case, it follows that to be 
stable a medium must satisfy 

®A%1 >° (28) 

The last two conditions are, in fact, equivalent. The condition for cy is trivially satisfied 
in a relativistic high temperature system, since the free energy density / is dominated 
by the term ~ — T 4 . Hence essentially only the last condition remains to be inspected. 

Let us note that there is an equivalent formulation for the second question. Indeed, 
the fluctuation matrix around the saddle point used in the evaluation of eq. ( |2.5| ) must 
be positive definite, from which it follows that 

^Wl«„>0. (2-9) 

To address the third question, one has to study non-homogeneous configurations in 
eq. (|2.5| ). To do this in full generality is, obviously, very difficult. However, when one 
restricts to spherically symmetric configurations (Q-balls), which is expected to be the 
relevant case, the problem becomes, to a good approximation, solvable. 

Finally, the fourth question is a non-equilibrium consideration, and the tools above 
do not directly apply. However, the equilibrium limit turns out to give quite strong 
constraints. As the results depend on the initial non-equilibrium state, we will not at- 
tempt a general analysis here, but rather concentrate on the cosmologically interesting 
case of Q-ball evaporation in the Early Universe. This will be done in Sec. ||]. 
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3 Q-balls at zero and finite temperatures 



It turns out that the first and second questions formulated in Sec. |2], concerning the 
existence and stability of Q-matter, have a simple and transparent answer: a Q-matter 
state (~ Bose-Einstein condensate) could in principle exist, but it is unstable and will 
decay into an inhomogeneous configuration. We will illustrate this phenomenon in 
some detail in Sec. f|. In this Section we consider the inhomogeneous final state, i.e., 
Q-balls, in general qualitative terms. Some more details in a specific model will again 
be given in Sec. |j. 

Let us start by reviewing the properties of Q-balls at zero temperature. Consider 
a generic field theory containing scalar fields <pi and fermions ipi with global U(l) 
charges qi. The scalar potential is V((p). Complications appearing with the inclusion of 
gauge fields were discussed in p7| and do not appear in the phenomenologically interest- 



ing case of SUSY theories, where flat directions are associated with SU(3) xSU(2) xU(l) 
singlets. Then, a spherically symmetric Q-ball solution has the form ||J-0 

4>i{r,t) = eycp(iqiUJt)(f)i{r). (3.1) 

The functions <pi{r) can be found by minimizing the functional 

E u = J d 3 x [£ |^-(r)| 2 + VL(0)], (3.2) 
i 

where 

VM = V(cj ) )-co 2 J2<lM<i>i- (3-3) 

i 

The frequency to is related to the total charge of the solution Q as 

Q = 2uJ d'x^^Ui- (3-4) 

i 

Consider now the finite temperature case and a large but finite volume V. The 
thermodynamics of the system is determined by the grand canonical partition function 
defined in eq. (|2.1|) . 

Given the Lagrangian of the system, which contains all the fields of the theory, one 
can derive a Euclidian path integral expression for Z^, 

Z^ = J V[A\\ fields] exp (- J^dr J d^xC^j , (3.5) 

where Ce is the ordinary Euclidean action with the replacement of the Euclidian time 
derivative [PB] 

d T -> d T - qiji. (3.6) 
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Note that for scalar fields, the charges come with the opposite signs in d T (j), d T <ft . As 
usual, bosonic fields are periodic, while fermionic fields are anti-periodic on the finite 
time interval < r < f3 = 1/T. 

To define if there are any non-trivial contributions to the partition function, poten- 
tially associated with Q-balls, one can look for the saddle points of the exponential in 
eq. (|3.5| ). The configurations that play the most important role at high temperatures 
are static, i.e., independent of Euclidian time. The static bosonic part of the action 
has the form 

s = pj ^x[^\^M 2 + v(4>)-^^44h], (3-7) 

j i 

precisely the one given by eqs. ( |3.2|) , ( |3.3| ) with the replacement u — > Thus, if the 
theory at zero temperature has Q-ball solutions with some values of u, these solutions 
are saddle points of the Euclidian finite temperature and finite density path integral 
at fi — u. The relationship between the charge and fi is now given by eq. ( |2.2|) , rather 
than eq. (pO|): besides Q-balls, at non-zero temperatures charge can be carried by 
particles outside it, as well. 

As is very well known from the study of phase transitions at high temperatures, the 
saddle points of the tree-level static bosonic action provide a good approximation only 
when T <C M, where M is the typical mass scale of the theory. At T > M, the quantum 
corrections are large, and must be resummed in some way. Thus one either has to use 
an approximation for the full effective action in eq. ( |2.5| ), or one has to use an 
effective theory where quantum corrections are small. The latter could be constructed 
for static but space dependent bosonic configurations, integrating out all fermions, all 
heavy (mass ^> T) bosons, and the non-zero Matsubara frequencies of light bosons. 
For static bosonic configurations, depending on space coordinates only weakly (we will 
see that this is indeed the case for sufficiently large Q-balls), the derivative expansion 
is most helpful: 

S cS = pJ d 3 x[j2Z j ((f>,T,^\d i( l )j (r)\ 2 + V eS ((l ) ,T,^ + ..], (3.8) 
j 

where V e ff((fi,T, fi) is basically the effective potential of the system at non-zero T and 
fj, (more precisely, V e s is the contribution to the effective potential from the degrees of 
freedom that have been integrated out). The effective potential as a function of \x has 
been computed in a number of theories and discussed in connection with the question 
of the influence of charge density on the symmetry behaviour in a number of papers 



28|-|30|1. If the particle masses in the background of the field <p are smaller than the 
temperature, then the high temperature expansion can be used for the construction of 
the effective action precisely along the lines of Refs. [|31|— 134|. If not (as is in fact for 



many of the considerations below), then an effective action incorporating both high 
and low temperature behaviour must be used. Now, we have at hands enough for 
a qualitative discussion and order of magnitude estimates of the properties of SUSY 
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Q-balls at high temperatures. For this aim we take the wave function normalizations 
in eq. fl3.7| ) to be unity. As a number of studies |35|, |31| show, higher order corrections 



in the wave function normalizations are not essential numerically. 

Let us take some field representing, in somewhat loose terms, a gauge-invariant 
combination of squarks and sleptons in the MSSM, and carrying baryon (and, perhaps, 
lepton) number (for simplicity, we take the charge q = 1). Specific examples can 
be found in HQ. Under the global baryon symmetry, the field is transformed as 
— ► exp(ia)<f). Assume that the effective potential is flat along this direction as shown 
in Fig. 1, 

V(0) -> M 4 (3.9) 

at —>■ oo. Then, at zero temperature the minimum of the energy at a fixed baryonic 
charge Q ^> 1 corresponds to a time dependent Q-ball solution which can in the limit 
of a large charge be written as 



(sin tor) /tor, r < Rq 
(r) ss exp(iwt) x { (3.10) 
0, r > R . 



Here 



R Q « (1/V2) M- X Q X I\ Q ^ (1/V2) MQ 1/4 , uR Q « tt, (3.11) 
and the mass of the solution is 



M Q ?3 (4ttv / 2/3) MQ 3/4 . (3.12) 

A Q-ball is stable (cannot decay into protons: M Q < m p Q) if Q > 10 15 (M/lTeV) 4 . 

Let us now see what happens at high temperatures. The scalar field couples to 
the other fields of the MSSM. In the background of this field a number of particles 
(e.g. gluons and gluinos) acquire masses m(0) ~ g<p, while others do not (e.g. if the 
Q-ball solution at T = is constructed from squarks only, leptons and sleptons remain 
massless since they do not have tree level interactions with squarks). Now, if g<f> > T, 
massive particles do not contribute to the finite temperature effective potential just 
because of the Boltzmann suppression exp(— m(0)/T), whereas the contribution of the 
massless particles is <fi independent and produces just an overall shift of the effective 
potential which we discard. At the same time, m(0) is small near = 0. Thus, the 
point = moves down as AV ~ —NT 4 , where the coefficient iV is related to the 
difference between the number of light degrees of freedom at small and large 0. So, the 
finite temperature effective potential has qualitatively the same form as the effective 
potential at zero temperature, with the change M 4 — > M(T) 4 ~ M 4 + NT A . 

This consideration allows to write down immediately the Q-ball contribution to the 
thermodynamic potential. Replacing Mq —>■ Fq, M — > M(T) in eq. ( |3.12| ) and using 
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that Q 1/4 ~ R Q M ~ M/u ~ M/fJL, we get 



FqM ~ M(T) , (3.13) 



with the charge in the Q-ball 

QM~(*f)'. (3.14) 

The size of the Q-ball is given by Rq ~ 7r//i. 

Now, the /i-dependent part of the thermodynamic potential has the form 

n(T, V, /i) = F Q (fi) - C^ 2 T 2 V, (3.15) 

where the coefficient C is related to the number of light carriers of the baryon charge 
(such as quarks). The total charge of the system follows from eq. ( p.2|) , 

Qtotal = Q Q (/i)+C>TV. (3.16) 

These two equations define equilibrium Q-ball properties at high temperatures. 

Let us assume first that the average charge density hq is nonzero and is fixed, while 
the volume of the system grows. Then eq. (|3.16|) admits two solutions. The first one is 
with a "large" chemical potential \i ~ hq/(CT 2 ). It corresponds to a state where all 
the charge is carried by particles in the plasma and Q-balls are not present. Another 
solution has a "small" chemical potential, 

/j, ~ M(T)(n Q V)~ 1/4: . (3.17) 

It corresponds to a system containing one Q-ball that carries almost all of the charge 
Qtotai = n Q V. 

Constructing the Legendre transform of eq. ( |3.15| ) according to eq. Q2.3|), it can now 



be seen that the latter solution has a smaller free energy at large volumes. Indeed, 
the Q-ball free energy scales with volume only as Fq ~ V 3 ^, whereas the free energy 
of the plasma phase grows as -F p i asma ~ V. Thus, we arrive at the conclusion that in 
theories with flat potentials and a fixed non-zero density of the charge Q, the ground 
state of the system always contains just one Q-ball at any temperature when V — > oo. 

Let us then take another limit, where the total charge of the system is fixed, and 
we increase its volume. In this case, the average density of the charge decreases as 
n Q ~ Q total/ V- Then, a state with a Q-ball is more favourable than a homogeneous 
distribution of the charge at 

V< Vtotal . (3.18) 
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This inequality allows, for example, an estimate of the temperature T evap at which an 
initially cold Q-ball, placed in a volume V, evaporates. It follows from eq. ( |3.18| ) that 
^vap ~ 1/- 1/3 Q 5/12 , provided T evap > M, and T evap ~ (VM)-^ 2 Q 5 ^ in the opposite 
case. 

All these results have a simple physical meaning. Imagine that we place a Q-ball at 
zero temperature into a box with the volume V. Now, let us heat the system gradually 
without adding any charge to it. Some charge from the Q-ball (5Q) will evaporate, 
and it will create a chemical potential in the surrounding plasma, /i p i asm a ~ $Q/(VT 2 ). 
The process of Q-ball evaporation will stop when the chemical potential associated 
with the Q-ball, \xq ~ M(T)(Q — 5Q)~ l l 4 , will be equal to /i p i asm a- If the inequality in 
eq. (|3.18|) is satisfied, then one finds that the solution for 5Q is 

jg ^ M(T)Vr> 

so that the Q-ball cannot evaporate completely, whereas if eq. ( |3.18|) is not true, then 
the equation yU p i asma = [Aq does not have physical solutions with Q — 5Q ^> 1, and 
therefore, the Q-ball disappears. 

As we see, the structure of the ground state is obviously volume dependent. Thus, 
the system we consider does not have a standard thermodynamic limit. The flatness 
of the potential is essential for these conclusions. If Fq ~ Q 1 rather than Fq ~ Q 3 ^, 
then the standard thermodynamical limit is well defined, and Q-balls cease to exist 
above some temperature which is entirely determined by the density of the charge and 
the parameters of the model. 



4 A model computation 

In order to illustrate in more specific terms the issues discussed above, we will in this 
Section consider in some detail the renormalizable model introduced in |J for studying 
super symmetric inflation. We first discuss the stability of Q-matter, and then that of 
Q-balls, adding numerical coefficients to the estimates in Sec. |[ 

The model is defined by a superpotential with two complex fields 0, x- 

W = l -f X <p 2 ~ Xm 2 . (4.1) 

The parameters f,m 2 can be assumed to be real. We denote m 2 = fm 2 . In the 
direction of the charged field, the potential of the model grows only logarithmically. 
Let us assume that there are no local symmetries and thus no gauge fields. 

We start by writing down the Minkowskian Lagrangian. Let us split <fi to real 
components, 4> — (0i +i<f>2)/V^, and denote the Weyl spinors corresponding to 4>, x by 
VV, il>x- Tnen 

C M = 5^ + + + - Vtfufax) - (4-2) 
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where 



Vfa,<h,x) = ^ 2 (^-0i) + i/ 2 (0i + ^) 2 + J/ 2 (0 1 + ^)W 2 + ! ^, (4.3) 



2 ^ u 16 y 1 ZJ 2 v 1 ^' 1 / 2 

Cp = -/x^ + /#x^+H-C (4.4) 

This action has a global U(l)-symmetry, 

4>^<P, X^ e ia X , 

^ -> e-^ 2 ^, ~> eW2 ^x- (4-5) 

The superpotential is not symmetric and therefore the scalar fields and the correspond- 
ing fermions transform differently. The charge corresponding to the symmetry is 



Q = Jd 3 x 



(4-6) 



The theory has also a discrete symmetry, — > —<p, ip^ — > —ip^, which however does not 
play any role in the following. 

At tree-level at zero temperature, the ground state of the system is at a broken min- 
imum: x = 0, 02 = 0, 4>1 = 4m 2 // 2 , V(4>, x) — 0. There supersymmetry is conserved, 
and the spectrum consists of four massive scalar degrees of freedom and one Dirac 
fermion, all with the same mass, y/2ih. On the other hand, for large values of 
the minimum of the potential is at <\>\ = <p 2 = and the tree-level potential does not 
depend on x : this is a flat direction. Along the flat direction, two scalar degrees of 
freedom and one Majorana fermion remain massless. The value of the potential along 
the flat direction is m 4 // 2 = m 4 . 



Consider then in eq. ( |2.1|) . According to the standard procedure pq] , the Euclid- 
ian path integral expression for with the charge Q as given in eq. ( |4.6|) , is 

Z» = V<f>VxV^V^ x exp (- J^dr J d^xC^j , (4.7) 

where, denoting now also x = {Xi + ^X2)/v / 2, the Euclidian Lagrangian is 
Cb = ^(^ 1 ) 2 + ^(^ 2 ) 2 + i(V0 1 ) 2 + ^(V0 2 ) 2 

+ \(drXl - WX2) 2 + \{drX2 + l^Xlf + l^Xl) 2 + \{VX2? 



+ V(<f ) ,x) + C x . (4. 
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4.0 



v i-ioo P (<l>i>2Ci>M)) 

(projection to the x,-axis) 



T/M=1.0 




Figure 1: The projection of the effective potential on the Xi-axis, at different tempera- 
tures and = 0. The temperature is given in units of the mass scale M = m — m(fh ). 
The thick lines indicate the regime where (0i) ^ 0: this symmetry is restored at 
T/m > 4.0. For Xi/T<3.5, the effective potential is the one of the 3d effective the- 
ory in eqs. ( |A.16| ), ( |A.17|) , and for xi/T> 2.5, the one of the 4d theory in eqs. ( |A.19|) , 
( |A.21| ). The scale of Xi is in both cases what appears in the 3d theory: xi = Xi d (/"r/4)- 



Here V(4>,x), £\ are fr° m e Q s - ( ^-3|) , 

We now wish to consider the functional integral for in the background of 0i, xi, 
i.e., to compute the effective potentials V^(0i,xi), Vq(4>i,Xi), defined in eqs. ( |2.6|) , 
( p.7|) . For X-, choosing to consider xi is no restriction, due to the U(l) symmetry. On 
the other hand, we have 02 = already at the tree-level at zero temperature, so 02 
remains at origin for T > 0. In the following, we denote 0i = 0i, xi = Xi- 

The details of the computation of the 1-loop effective potential are discussed in 
Appendix A. We consider here the numerical results and their main features. To have 
a weakly coupled theory, we take / = 0.5. 



4.1 The stability of Q-matter 

In Fig. m the 1-loop effective potential is shown at fi — 0. It is seen that for T/m ^,4, 
the symmetry is restored (0i = 0) for any Xi- However, irrespective of the temperature, 
the effective potential always displays the characteristic form allowing Q-balls: at small 
Xi, there is a thermal or vacuum mass term, but at large xii the potential flattens off 
and grows eventually only logarithmically. 

Consider then the case of a finite charge density. In Figs. 0, |3|, the potentials V^(xi) 
and its Legendre transform V q (xi) are shown at T/rh = 5.0,2.0, respectively. The 
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T/M=5.0 T/M=5.0 




5C/T 5C/T 

Figure 2: Left: The effective potential V^(xi) at T/m = 5.0, for different values of 
\ijT . Right: The effective potential V q (xi) at T/m = 5.0, for different values of q/T 3 . 
Each extremum of V^(xi) corresponds to an extremum of V q (xi), as well, as is shown 
by the symbols. For q/T 3 > (q/T 3 ) c ~ 0.065, there is a charged condensate: (xi) > 0. 



structure is simpler at T/fho = 5.0, where the 0i-symmetry is restored everywhere and 
the mass scale m does not have much effect. Nevertheless, both cases show the same 
general behaviour. There is a critical charge density q c > 0. For q < q c , the minimum 
of V q {xi) is at xi = and the charge is hence carried by particle excitations in a plasma 
phase. However, for q > q c , a phenomenon analogous with Bose-Einstein condensation 
takes place, and xi > 0. 

This pattern can be understood already from the tree-level potential of the effective 
3d theory in eq. ( |A.16| ). Consider first the temperatures T > « 2m/ f ~ 4m. Then 



ml > according to eq. (|A.15|) , and (pi = for all Xi- The effective potential V^(xi) is 



V,(xi) ^ ^st. - |/, 2 T 2 + \ (V + £pj x\. (4.9) 

It will be useful to consider a more general expression for the latter term in the coef- 
ficient of x\i so let us replace f 2 T 2 /8 — > M 2 . Then the situation looks precisely like 
relativistic Bose-Einstein condensation in a free theory. 

To proceed, let us make a Legendre transformation into V^(xi) = ^t(xi) where 
q = -dftVftixi)-- 

v «~\m&T7i + ¥ 12 ^ (4 - 10) 
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T/M=2.0 T/M=2.0 

-0.40 i . 1 . 1 . 1 i . 1 . 




Figure 3: As Fig. |^, but for T/rfiQ = 2.0. The thick lines indicate the regime where 
the symmetry is broken, ^ 0. Due to the symmetry breaking, the potentials have 
more structure than in Fig. |^. The critical charge density is here (q/T 3 ) c m 0.25. Note 
that the (unique) extremum of V^(xi) at Xi > is always a maximum. 



What then remains is to minimize V q (xi) with respect to xii f° r different q, T. The 
main results are as follows. The global minimum is at xi — if M > (12/5) (q/T 2 ). 
Hence, the critical charge density is 



q c = j^MT 2 . (4.11) 



If q < q c , then 



Xi = 0; H=~l /(r,ff) = /(T,0) + 2i; > 0. (4.12) 

If q > q c , on the other hand, then 



Xi 



q -J±- fi = M- f(T,q) = f(T,0) + ^(2q-q c y, ^(Xi) = 0. (4.13) 



Replacing now M 2 — > f 2 T 2 /8, we see that for any given T > Tf ~ 2m/ f, there is a 
critical charge density 

±) (4 14) 

X'Jc 24yT 1 } 
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0.30 




o.oo L-^ — 1 — ■ — 1 — ■ — 1 — ■ — 1 — ' — 

0.0 0.2 0.4 0.6 0.8 1.0 



Figure 4: The free energy density as a function of the charge density for T/rho = 
2.0,5.0. For T/m >4 the symmetry is always restored (0i = 0) and the results scale 
with T; the critical charge density below which xi = 0, is (q/T 3 ) c « 0.065. For 
temperatures T/rho ^,4, on the other hand, the symmetry is broken (0i > 0) at small 
Xi, and there is more structure. 



below which there is no condensate, Xi = 0- Note that for T > Tjt, the leading 
order results in eqs. ( |4.12|) , (|4.13| ) do not depend on fh at all, and all the dimensionful 
quantities scale with T. Numerically, (q/T 3 ) c ~ 0.07, in accordance with Fig. 0. 

For T<T^, on the other hand, m\ < and <p 2 > for small xi- The potential 
in eq. ( A. 16 ) must thus first be minimized with respect to <fi\. After that has been 



done, the potential V^(xi) is again of the form in eq. ( |L9p for small xi, but now with 
M 2 = 2m 2 — 3f 2 T 2 /8. Thus the mass scale m appears in the results. For example, 



t\ 5./a*> w t (415) 



T 3 J C 12 V T 2 

which gives (q/T 3 ) c ~ 0.27 for T/m = 2.0, in agreement with Fig. |3|. 

To summarize, based on free energy considerations alone, one would say that for any 
temperature, there is a charge density q c such that for q > q c , one has a Bose-Einstein 
condensate, or a "Q-matter" state. For T — > 0, q c — >• according to eq. ( |4.15| ), and the 
plasma phase does not exist at all. 

Let us then consider whether such a Q-matter condensate would be stable. Apply- 
ing the stability conditions in eqs. ( p.8|) , (|2.9| ), one sees that in the tree level case, 
the plasma phase [eq. (|4.12|) ] is stable while the situation in the condensate phase 
[eq. ( [4.13|) 1 is marginal: d 2 f(T,q) = d^V^Xi) = 0- Hence the system is stabilized or 



destabilized through interactions. Adding an interaction of the type \(x*x) 
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in fact stabilizes the condensate. On the other hand, it is easy to see that in the present 
case, the interactions destabilize the condensate! 

The fact that the condensate in the present system is unstable, can be observed from 
Figs. |3|. The extremum corresponding to the condensate is always a maximum of 
V M (xi), and thus the inequality in eq. ( |2.9| ) is violated. Equivalently, one can observe 
from Fig. |] that the free energy density f(T, q) has a negative curvature at q > q c , which 
violates the stability condition in eq. ( |2.8| ). The negative curvature can be understood 
also analytically: at large charges, corresponding to large xi's, 

vAx.) fx £ T v,(o) + 5(i- 

■2 T , 2 T 2 , 1 m 2 

+ ^T + -W-2^ + ^ fn ^W (416) 




from which it follows that 

f{T,q) /(/■.()) 4 ^ 




+ -^r + ^™H _ * . (4.17) 



Hence dgf(T,q) < 0. We conclude that a (meta)stable Q-matter state does not exist 
in the present system. 

Consider then the plasma phase. That the plasma phase at q < q c satisfies the 
stability conditions in eqs. Q2.8|) , (|2.9[) , only proves that such a state is metastable 
against small perturbations. The global minimum may still be elsewhere. Consider, in 
particular, an arrangement where all the charge is moved into a subvolume V = rV, 
with r < 1 (this of course means that there is no chemical equilibrium, but this is 
not essential for the argument). Let 5f(T,q) = f(T,q) — f(T, 0). Then the original 
free energy density excess with respect to an empty space is 5f(T, q), while the final is 
5f(T,q/r)r. Now, for any initial q, 5f(T,q/r)r — > as r — > 0, since 5f(T,q/r) grows 
only logarithmically at small r, see eq. ( [4.171 ). Hence, even if f(q/r)r grows at r < 1 
and the original state is thus metastable, f(q/r)r finally starts to decrease for small 
enough r and it is in principle favourable to clump all the charge together in one place. 
Of course, surface effects start to be important at very small r, so to be more precise, 
one has to inspect the possible final configurations in more detail. We thus turn to 
Q-balls. 

4.2 The stability of Q-balls 

We start by considering the limit of large temperatures, T 3> rh, and small chemical 
potentials, |i <C T. Moreover, let us subtract the contribution of the homogeneous 
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symmetric plasma phase, where xi = 0- Then it follows from eq. (|A.21 ) that 

1 



U (MJT) - -U2M X /T) 



(4.18) 



where M x = fxi/ V% and Jo is from eq. ( |A.22j ). Neglecting the change in the derivative 
term, the extra free energy of any configuration Xi( r ) is then given by 



d 3 r 



Mxi) 2 + v(xi) 



(4.19) 



The properties of the saddle points of eq. ( |4.19| ) are discussed in Appendix B. It 
follows from eqs. flBJf), @3g ) that in the present model (a v = (vr 2 /24)T 4 ), the Q-ball 
free energy, charge and radius at T>m are given by 



F Q « 4.74TQ 3/4 , Q « 160 



R 



7T / /i 



T VT 



(4.20) 



We now consider the stability of such Q-balls. The other metastable alternative is 
the homogeneous plasma phase, so we have to compare with f(T, q) in eq. (|4.12j ). Let 
us write F Q = aTQ 3/4 , Q = &(/i/T)" 4 . 

Consider first a situation of thermal but not chemical equilibrium, in which all the 
charge is in the Q-ball. Then, the Q-ball does not decay into a homogeneous plasma 
phase, if F p i asma = V/ p i asma (T, q) > F Q , or 
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- 4 Q 5 / 4 > VT 3 . 



(4.2i; 



This is the analogue of eq. ( 3.18| ), and tells that in a fixed volume, Q-balls will evaporate 
at large enough temperatures. 

On the other hand, note that what appears on the right-hand side of eq. Q4.21 ) is 
a comoving volume in units of temperature, which is a constant in the cosmological 
context. Thus, in this case, if Q-balls are thermally stable at some high temperature, 
they remain stable at least as long as T > m. The thermal stability of a distribution 
of Q-balls can be seen by replacing V by their inverse number density rig 1 . 

Consider then a thermodynamical equilibrium situation where there is chemical equi- 
librium, as well. Some of the charge will now be outside the Q-ball, since there is a 
non-vanishing chemical potential. The stability condition becomes 



6 Q 2 



(4.22) 



where the Q-ball charge Q' is to be solved from an equality following from chemical 
equilibrium: 



pla 

T 



12 Q-Q' 
5 VT 3 



Q-ball 



b 



-1/4 



(4.23) 
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However, it can easily be seen that when eq. ( j4.21|) is strongly enough satisfied, then 
the fact that there is chemical equilibrium does not change the result. Indeed, it can 
be seen from eq. ( }4.23|) that most of the charge resides in Q-balls in this limit, 



2«1 

Q 



~2a 



VT 3 



(6/5)0-^5/4 



(4.24) 



and Q-balls carry most of the free energy associated with the non-zero net charge, 



plasma 



&1/2 
4^2 



VT 3 



<C 1. 



(4.25) 



h q L (6/5)a" 1 g 5 /4_ 

Moreover, provided that the Q-ball charge is large, Q-balls are small in the sense that 
most of the space is in the plasma phase: 

V Q (A/3)tt(R q T) 3 10vr 4 a \ (6/5)^^^ 



plasma 



if 



VT 3 



Q > 130 



96 3 / 4 



VT 3 



QV 2 



<1, 



(6/5)a~ 1 Q 5 / 4 
VT 3 



(4.26) 



(4.27) 



Even though chemical equilibrium does not essentially change the constraint in 
eq. ( |4 . 2 1| ) , there is an important implication following from eq. ( ^4.23| ). Noting that 
the entropy density is s = dp/dT = n 2 T 3 /3, one gets a simple relation between the 
charge density in the plasma, and the charge residing in a Q-ball: 

q\ 5 5i> 1//4 

A _Q 

plasma 



4tt 2 T 4tt 2 



-1/4 



0.45(Q') 



1/4 



(4.28) 



In the cosmo logical context, this corresponds to a relation between the baryon asym- 
metry in the plasma, and the baryon number residing in Q-balls. 

Finally, consider what happens at lower temperatures. The Q-ball free energy will 
decrease for a while as ~ T, until at T<rh/f, the symmetry breaking starts to play 
a role: eventually T gets replaced by m// 1 / 2 in eq. (|4.20| ). At the same time, the 
Q-part of the free energy of the plasma phase also decreases, initially as -F p i asma ~ 
(6/5) (Q 2 /VT 3 )T ~ T. This behaviour is valid up to T ~ m/ir; at smaller temperatures 
one gets the standard non-relativistic ideal gas behaviour -F plasma ~ QT[\nQ/(VT 3 ) + 
(3/2) ln(T/m)] ~ T. Hence the plasma phase free energy decreases faster than Fq 
and a thermalized Q-ball in a comoving volume becomes more unstable at smaller 
temperatures. 

The considerations above concern an equilibrium situation. There are also other 
relevant questions, such as, assuming that one is trapped in a metastable homogeneous 
plasma phase and Q-balls would be stable, what is the transition rate? Some progress 
in addressing this question has been reported in f3"6fl , but we will not consider it here. 
In the next Section we consider the opposite limit: taking a Q-ball which is unstable, 
at which rate does it decay into the plasma phase? 
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5 Baryon to cold dark matter ratio 



Consider Q-balls in the cosmological context. In || it was shown that Q-balls with 
large charges can be produced in the Early Universe. The general idea is that in 
the AD scenario for baryogenesis, an uniform condensate of the scalar field, carrying 
baryon number, is formed as a result of inflation, CP-violation and baryon number 



non-conservation fll2| . This AD condensate is nothing but Q-matter, described by 
an uniform scalar field with a time dependent phase. As we have seen, Q-matter is 
unstable against Q-ball formation; an analysis of the typical instability scales shows 
that charges of the order of magnitude Q > 10 20 can be easily formed (and do not 



violate experimental constraints p7|). In || it was suggested that these stable Q-balls 
can play the role of the cold dark matter. The ordinary, baryonic matter appears then 
as a result of Q-ball evaporation. If true, baryonic and cold dark matter in the Universe 
come from the same source. Moreover, dark matter is in fact baryonic in this case, but 
baryons exist in the form of SUSY Q-balls which do not participate in nucleosynthesis. 
In this scenario the ratio of baryonic to cold dark matter can be computed and the aim 
of this Section is to make such an estimate. 

We assume that after the instability has taken place, the Universe is filled with 
baryonic Q-balls with charge?] Q and number density uq, surrounded by a hot plasma 
at temperature T. We will also assume that only a negligible part of the baryonic 
charge initially exists in the plasma. Some arguments if favour of this assumption were 



put forward in [25], but its complete check would require the solution of the non-linear 
problem of Q-ball formation. In the following, we completely neglect the process of 
merging of two Q-balls, since their concentration is so small that the probability of 
their collision is absolutely negligible, as can be checked a posteriori. 

As a starting point, let us consider what the thermodynamical equilibrium state 
corresponding to these initial conditions would be. Then we will argue that in fact 
only thermal equilibrium can be established in practice, while chemical equilibrium 
with respect to baryon charge is not reached with these initial conditions. 

If there is full thermodynamical equilibrium, then it follows from eq. ( |4.23| ) that 
already at high temperatures, the baryon to photon ratio in the plasma phase would 
be 

v = — - Q- 1/A - (5-1) 

In full thermal equilibrium Q-balls become more unstable at lower temperatures, see 
Sec. P~2, so that with the experimental number rj ~ 10~ 10 , eq. (|5.1|) would imply the 



lower bound Q > 10 . On the other hand, if Q-balls account for a fraction f2 DM of the 



3 As is shown in ||, the instability of an AD condensate occurs in a narrow range of scales, and 
thus the Q-ball charge distribution is likely to be narrow, as well. 
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dark matter, then their present number density is tSq = Pq^/Mq, where 

p§ « 10- 5 fi DM m p cm" 3 , (5.2) 
and m p is the proton mass. It follows that 

^T^3xl0 8 ^^, (5.3) 

Wlp i 'DM 

where g* is the number of massless degrees of freedom. If, for instance, g* ~ 200, Q~dm ~ 
0.4, Q ~ 10 40 and M ~ (1 . . . 10) TeV, thenf] n^T 3 ~ 10 39 . . . 10 40 . Comparing now 
with the LHS of eq. (|4.21j) , it is seen that such a density of Q-balls would indeed be 



thermodynamically stable against evaporation into the plasma phase. However, it is 
not clear how Q-balls with charges as large as 10 40 could be produced, and whether they 
would indeed reach thermodynamical equilibrium in the cosmological environment. 

Now we are going to take into account the expansion of the Universe. First of 
all, the intrinsic temperature of Q-balls may be different from the temperature of the 
surrounding plasma, if the system is not fully thermalized. Then Q-balls are heated up 
because of collisions with particles from the plasma. The energy transfer to the scalar 
particles in the condensate from which the Q-ball is built, has been estimated in || 
(we assume that all coupling constants are of order unity), 

d j- t ~ R 2 Q T 3 an,6R6E ~ T 2 (^) 3 Q 1/4 - (5.4) 

Here RqT 3 counts the initial flux of particles in the plasma on the Q-ball, 5E ~ uo 
gives the energy transfer to a scalar particle in the condensate, 5R ~ T/M 2 gives the 
thickness of the Q-ball layer in which plasma particles can penetrate (due to ~ 
^I^T), n 4> ~ 4> 2uJ ~ T 2 MQ~ 1 / 4 is the charged scalar particle number density in this 
layer [eq. (|3.4| )1, a ~ l/(uT) is the typical cross-sectionF] of the interaction of an 
energetic particle from the plasma with a condensate scalar with a small energy ~ u, 
and we assumed T 3> M. The estimate in eq. ( |5.4|) contains only the effect of particles 
interacting strongly (at tree-level) with the condensate, so that it should represent a 
lower bound of the thermalization rate. 

The energy transfer is most efficient at high temperatures. Requiring that the change 
of Q-ball free energy is of the order of its equilibrium value F ~ TQ 3//4 , we find that 
a zero-temperature Q-ball can be thermalized if the initial temperature of the plasma 
satisfies the requirement 

/ M \ 1/2 

T in >M( — ) Q'l\ (5.5) 



4 To illustrate this number density in cosmological units: ng 1 ' 3 ~ 10 2 R™, where is the 
horizon radius at the electroweak scale T ~ 100 GeV. 

D In m this cross-section was incorrectly taken to be a ~ 1/T 2 . Also, the estimate of the transmitted 
energy SE contains a misprint: the power 2 must be omitted from the corresponding expression. 
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where M = (0.30/ v /^)m P1 ~ 3 x 10 17 GeV. For typical values of M ~ 1 - 10 TeV 
and Q ~ 10 20 — 10 30 , Ti n > 1 . . . 10 4 GeV, which is true for most inflationary models. 

Next, we would like to determine the amount of baryons which evaporated from 
thermalized Q-balls. From the discussion in Sec. |3|, we know that we can associate 
with Q-balls a chemical potential \x ~ M(T)Q~ 1/4 . In the case that the surrounding 
plasma has the same chemical potential, Q-balls are in chemical equilibrium and their 
charge does not change. In other words, the rate of their evaporation coincides with 
the rate of charge accretion on the Q-ball. The upper limit on the accretion can be 
easily computed: the rate of accretion cannot be larger than the total baryonic flux 
through the Q-ball surface. So, 

-fa = -D((J.Q -/^plasmajT AttR q ~ -D , (5.6) 

where the coefficient D < 1 shows the deviation of the real evaporation rate from the 
upper bound. At T > (m^ is the mass of squarks, suppressed by some coupling 
constants with respect to M), the evaporation rate must be parametrically the same 
as the upper bound, since the scalar particles building up the Q-ball appear with 
a large number density ~ T 3 also in the plasma phased Hence, the rate of Q-ball 
evaporation in a plasma with a vanishing chemical potential (/i p i a sma = 0) at T > 
is given by eq. (|5.6|) with D ~ 1. At T < m^, the charge in the plasma sits mainly in 
light fermions; the probability of squark emission from a Q-ball is suppressed by the 
Boltzmann exponent exp(— m^/T). Therefore, the emission of heavy squarks from the 
Q-ball surface can be neglected. The main process which is allowed is squark-squark 
transition into two quarks, 00 — > qq. It may occur through gluino exchange with a 
cross-section of the order of a ~ (3 /mi (see also a similar discussion in H38H), where 
(3 is some combination of the coupling constants. Thus, we expect D ~ f3T 2 /m 2 ^ at 
T < m«£. 

By integrating eq. ( |5.6| ) (we assume a radiation dominated Universe) we find that 
Q-ball evaporation is negligible at T < (D is essentially replaced by (3 in eq. (|5.7|)). 
The charge Qt of a Q-ball at T ~ is related to its initial charge Q as 




^n 1 - 1 ^" 1 ^ 1 ' (57) 

Most of the charge is evaporated at low temperatures near T ~ m^. 

From eq. ( f).7\ j we conclude that Q-balls survive evaporation if their initial charge is 



6 In case some of the Lagrangian squark mass parameters are small or negative, this regime extends 
down to the critical temperature of the electroweak phase transition. 
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If, for example, M ~ 1 TeV, then all Q-balls with Qq > 10 20 stay till the present time. 
It is interesting to note that these values of Q are naturally produced from the decay 
of an AD condensate. 

Now, if eq. Q5.8| ) is satisfied, then the amount of charge evaporated from an individual 
Q-ball is of the order of 

AQ~D*£H— )(&'*■ (5-9) 

Thus, the number density of baryons ns in the plasma is related to the number density 
of Q-balls tiq as 

n B ~ n Q AQ. (5.10) 

Using then rig from eq. ( [5.3D together with the entropy density s pd (2/45)g*7r 2 T 3 , and 
assuming that Q-balls give all the cold dark matter in the Universe with Qbm ~ 0.4, 
the baryon to photon ratio is computable: 

n B 1n _ 8 M m P , ( M \ fr ,i\ 

V = 10 : ^T72 ln ( )• (5-11) 



n 



7 



The correct ratio i] ~ 10" 10 appears quite naturally with the choice of, e.g., M ~ 1 — 10 
TeV and a corresponding initial charge Qq ~ 10 28 — 10 22 . These values are quite 
plausible from the mechanism of Q-ball formation. 

Since Q-balls interact very rarely, they should not take part in light element nu- 
cleosynthesis. From the point of view of structure formation, they are completely 
equivalent to weakly interacting cold dark matter. They thus behave in a somewhat 
similar way as quark nuggets [ [39(1 . 



6 Conclusions 

Supersymmetric non-topological solitons in theories with flat directions of the effective 
potential have a number of interesting properties at high temperatures. Provided that 
their charge is large enough, Q-balls can exist at any temperature. Unlike for "classic" 
Q-balls with Mq ~ Q 1 , the ground state of a system with a finite density of the 
charge is not translationally invariant and contains Q-balls. The state of Q-matter 
(Bose-Einstein condensate) appears to be always unstable and it decays into Q-balls. 

These properties have interesting implications for cosmology. Depending on the 
mechanism of Q-ball formation, thermal and chemical equilibrium may or may not be 
reached. If it is reached, an equilibrium distribution of Q-balls with charges Q > 10 40 
could account for both the net baryon asymmetry needed for nucleosynthesis, and for 
the cold dark matter. In this picture, the dark matter is in fact baryonic, with the 
baryon charge distributed between ordinary baryons with n B /n 1 ~ 10~ 10 , and squarks 
packed inside dark matter Q-balls. 
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On the other hand, for charges Q ~ 10 22 . . . 10 28 which are more likely from the point 
of view of Q-ball formation, only thermal equilibrium is reached. But even then, it turns 
out that the right orders of magnitude for the baryon asymmetry in ordinary baryons 
and for the dark matter in Q-balls, can be reached provided that supersymmetry is 
broken at a low energy scale M ~ 1 ... 10 TeV. Thus Q-balls represent an appealing 
candidate for the cold dark matter, offering a new picture in which the origin of ordinary 
baryons and the dark matter is the same. 
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Appendix A: A model effective potential 

In this appendix, we derive the 1-loop effective potential of the theory considered in 
Sec. f|for given [i, T, both at small and large values of the fields. These two cases have 
to be treated separately, as the theory requires a resummation for small values of the 
fields. 

At 1-loop level, one has to consider renormalization. To be specific, the running of 
the fields and parameters is at 1-loop level determined by 



Tl-pf = t^-J\ JL^ = -^- 2 fm\ (A.l) 

d 2 „ . d 



= -is?^* ~^ x ' x = -mr> fx " x - (A - 2) 

Here 0*0, x*X denote the renormalized fields squared appearing in the Lagrangian, 
and not any composite operators. As a result, the parameters and fields at scale ~p are 

2/ 2 f 2 

« 0*0Ou o )(-^) ^ , X *x(-p)~X*x(Tio)(^-) , (A.3) 

where A = ~p exp[167r 2 /5/ 2 (/Io)]- 

The theory can then be defined in terms of two parameters, the values of m 2 (/I ), 
f(po), at some value JZ of the MS scale parameter. We choose /J = m , where 
m = m{mo). Then all the dimensionful quantities, such as the temperature T, are 
measured in units of m , and the results only depend on /(mo). To have a weakly 
coupled theory up to very high scales, we arbitrarily choose /(ttiq) — 0.5. 
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The most convenient way to implement the resummation needed at small values of 
the fields, is to construct an effective 3d theory describing the thermodynamics - 



34j. On the other hand, for large values of Xi no resummation is needed and the 
effective 3d theory is not valid any more. Then one can use the standard 4d finite 
temperature 1-loop potential. Let us discuss the 3d theory first. 



A.l The 3d effective theory and V((f),x) a ^ small 0, \ 



Let us denote 



Ma 



1 

71 J 



n, 



V2 



fXi- 



(A.5) 



The effective theory constructed is valid for M^, M x -C irT, and implements the re- 
summation needed in the regime x < fT. 

The derivation of the effective theory proceeds by computing the n-point functions in 
the original theory and in the effective theory, and by matching the parameters of the 
effective theory such that the original results are reproduced. In the present context, 
one has to consider the 2-point functions of the bosonic fields at non-zero momenta to 
get the field normalizations, and the 2- and 4-point functions at zero momenta. The 
field normalizations arise from fermionic loops alone. The zero-momentum correlators 
arise both from bosonic and fermionic loops, and the most convenient way of deriving 
them is the effective potential, expanded in terms of masses to the required order. We 
will also expand in /z, assuming /i <C ttT, although this would not be necessary. We 
then assume that parametrically f(p, fx, ^ ~ fT and work at 1-loop level to order f 4 . 
Thus we need the first five terms (up to B A ) from the expansion 



- In det{A + B) = - In det A 
2 2 



E 

n=l 



-1 



in+1 



2n 



-Tr (A B) n . 



(A.6) 



In the bosonic case, the matrices are 



A' 



B 



P 



diag(l, 1, 1, 1) 



/ m 2 + \Ml + M\ 



V 









-m 2 + |M| + Ml 
2M,,M X 






2M 4> M X 
V + Ml 



2fip, 



o 






-2/ipo 



(A.7) 



where p 2 = Pq + p 2 and po is a bosonic Matsubara frequency. In the fermionic case, 




B 



■f 7 o7 5 + M x M 

Ma, f 7 o75 
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The fermionic matrices have been written in the basis of the Major ana fermions 



^6 



(A.9) 



where we used the standard notation of Ref. ||40|| . There is an overall minus sign in 
eq. ( |A.6| ) for fermions, and only even powers of B contribute due to the momentum 
integration. 

The basic integrals appearing in the evaluation of the traces are 



'vbV 



12' 



p b {p 2 ) 2 167T 2 



7 + tw, 



C(3) 



'pfP 



24' 



(p 2 ) 2 



16tt 2 



> Pb {p 2 Y 128tt 4 T 2 ' 
1 r-\\ ^ 1 7 C( 3 ) 

7 + L/(/i) J' tW"i^¥ : 



(A.10) 
(A.ll) 



where Pb,Pf are the bosonic and fermionic Matsubara momenta, a prime means that 
the zero Matsubara mode is omitted, and 



L b (jl) = 2 In JL, Lfftl) = 2 In J±, Ju T = Ane^T. 



(A.12) 



The other integrals appearing, of the form % {Po) n / {p 2 ) m , can be derived from those in 
eq. (|A.11|) by taking derivatives with respect to T. 

As a result, the 3d finite temperature effective theory describing the thermodynamics 
of the original 4d theory is defined by the action 



S: 



3d 



V 

T 



+ J d 3 x 



tt 2 n 7 

. 7 Lt a --^ 2 T 2 + —-u 4 + 
12 24 P 192tt 2 ^ 



/ 2 (Mt) 



i 2 j.2 i 2 j.2 i 2 * 

+-m 1 1 + -m 2 2 + m x x X 
+ ^A (0 2 + 2 ) 2 + iA m (0 2 + 2 )x*x + A X ( X *X) 2 



(A.13) 



Here the fields are related to the 4d fields by 

1 



3d 



1 

T 



id 



^t). (XX)3d 



The parameters are given by 



(A.14) 
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(2- 
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16tt 2 
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24 



m 2 



(2- 4/3 ^) - 
P N 



/V f 2 T 2 
16tt 2 4 : 



16tt 2 



It/ 2 (2 8 / 5 /I r ) , X m = T/ 2 (2 2 / 5 7I T ) 



A, 



f 4 

8tt 2 



(A.15) 



Note, in particular, that around the two "symmetry restoring" temperatures [T^ = 
2m/ f where m\ ~ 0, and T c x = 2\/2/V/ where m 2 ~ 0], one has nT ^> {m, fi} for a 
small coupling /, so that the construction of the effective theory is well convergent. 

To go further with resummed perturbation theory, one can compute the effective 
potential in the theory defined by the action in eq. ( |A.13|) . The general pattern can be 
seen already from the tree-level potential, 

Vtrw, = const. + ^m 2 ^ + ^m\x\ + ^\/><f>i + ^m^lxl + jKxt- (A.16) 

This tree-level potential of the 3d theory corresponds to the dominant thermal screening 
contributions in the 4d 1-loop potential. The 3d 1-loop correction is 

1 



loop 



E 



12n . f ," J ' 

*=<P2,X2,± 



(A.17) 



where 



m A 



rn 



± 



ml + 



\ + ^ m xi, rn 



2 

X2 



m li + m L ± \/( m li - mlJ 2 + 4m. 



rn 



4>x 



m\ + 3X^1 + -X m xl; 
A TO 0iXi- 



m 



Xl 



(A.18) 



The 2-loop corrections can also be readily computed, but they do not affect our con- 
clusions. 



A. 2 The effective potential V((f), x) a t large % 

At large field values, M x > ttT, the 3d effective theory breaks down. Then one needs the 
full finite temperature 1-loop effective potential, but on the other hand resummation is 
not needed. Moreover, the ^-symmetry is restored, = 0. We assume also that the 
chemical potential is small, fi <C M x , which is well justified in the region considered. 
The tree-level potential for = is 

v *~ = TM ~ V*- (A ' 19) 
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The mass spectrum entering the 1-loop potential consists of two massless scalar degrees 
of freedom Xi,X2 and one massless Majorana spinor ip x , together with the massive 
scalar particles 0i, 02 and the Majorana fermion ip^. The masses are 



m ^ = M xi m 0i = ~ m2 + m 



(A.20) 



To order (n/M x ) 2 , the renormalized 1-loop contribution is then 



V x - 



loop 
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1=1,2 

(-2) 



8/ 90 48 P 
/I 2 3^ 



64tt 2 

,4 



In 



m 



^4> 



m l 2 
7 2 



+ - + 



2^2 



In 



64^ I m 2 ' 2 



+ ^ + 



ji4 



2tt 2 V8 



-J {2m f ,/T) - J {m^JT) 



In 



16tt 2 " m 2 ^ 



+ 



T 2 



2tt 2 V2 



-M2m^jT) - I Q {m^/T) 



(A.21) 



The integrals appearing here are 

J (y) = J^dxxHn (l -e-v 7 ^ 

ib(y) = 



dx 



\Jx 2 + y 2 e \/x 2 +y 2 



(A.22) 
(A.23) 



In the regime M x < T, the convergence of the 1-loop potential can be improved by a 
resummation: 



vi 



loop 



Vx —loop 



i=l,2 



T 
127T 



(A.24) 



where m\. = m\. + / 2 T 2 /4. 



Appendix B: Q-ball properties 



In this appendix, we briefly review the properties of Q-balls in theories with flat direc- 
tions. The high temperature limit of Q-balls in the theory considered in Sec. ^ follows 
as a special case. 

Consider the action for the length xx of the complex U(l) scalar field x (the extremum 
solutions considered do not depend on the phase of x), 



d 3 r 



l(diXi) 2 -\» 2 X 2 x + V(xx) 



(b.i; 
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Here it has been assumed that the derivative part of the effective action remains the 
same as at tree-level and only the potential changes. The potential is assumed to 
be normalized such that V(0) = 0. Let us scale eq. ([B.l ) into a dimensionless form 
by V(xi) = a v g(x), Xi = a x X, V = V^ r = a R r } F = a F F,Q = a Q Q, with = 
O'R — ®x/ y/ Q'V i ®f — \J &q — ■ The scaling of V(jfx) is chosen so 

that g(0) = 0, g(oo) = 1. It then follows that the spherically symmetric extrema of the 
action in eq. ( p.l| ), required by eq. fl2.5|) , satisfy 

^ + ~ = -^ 2 x + g\x); f(0) = 0; x(oo) = 0. (B.2) 



The main properties of the solution of eq. ( |B.2| ) for small fx are familiar from the 
zero temperature context, and do not depend at all on the functional form of g(x) at 
small X- To see this, assume that the solution scales so that = Cji~ a . Then 

r BF 

Q = l~lJ d 3 rx 2 = = Caft- a -\ (B.3) 

and 

F Q = F pL + ^Q={l + a~ l )(aC)^Q^. (B.4) 

The parameters a, C can now be easily derived. At large x (small f), the potential 
g(x) equals unity, and the equations of motion can be solved analytically: 

_ sin/if 

X{r) = Xo— — ■ (B.5) 
fir 

The scale of xo is determined by when — (l/2)/i 2 x 2 compensates for the unity in g(x)- 

It follows that xo ~ P> ■ The radius of the Q-ball is determined by x(R) ~ 1, i.e. 
R ~ irji^ 1 . Thus 

Q = fi[ d 3 rx 2 ~ /IR 3 Xo ~ r\ (B.6) 



and it follows from eq. ( |B.3| ) that a = 3. To determine the constant C, note that 
g(n/p) — > 6{k) with 8(0) = 0, when pL — > 0. Eq. ( |B.2| ) can then be integrated from 
R — e to R + e, and one finds xo — V2ft/jj" One then gets from eqs. (|B.3|) , (B.5| ) that 
C = 4vr 4 /3. 

We can now go back to the original dimensionful units, to find that 

R = (1/V2>y 1/4 Q 1/4 , xr x = a^ /4 Q 1/4 , Rli = vr, (B.7) 

F Q = (47rv / 2/3)a^ /4 Q 3 / 4 , Q = ^a V fi~\ (B.8) 

These are in accordance with eqs. ( |3.11|) , (|3.12|) , where ay = M 4 . In the model of 
Sec. (| at very high temperatures, ay = (7r 2 /24)T 4 . 
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Figure 5: The free energy and the charge Q as a function of the scaled chemical 
potential /2 (denoted by F', Q', /z' in this figure) for the Q-ball solution, in the limit 
of high temperatures and a small chemical potential. The extrapolated constant at 
Ji —* is C ~ 129.87(1), in accordance with the analytic result C = 47r 4 /3. The slopes 
are such that F M = C/r 3 (l - 2.1/2), Q = 3C/r 4 (l - 1.4/2). 



Finally, note that in contrast to the leading terms, the first corrections to the asymp- 
totic behaviour in eqs. ( |B.7| ), ( |B.8|) depend on the detailed form of the potential g(x)- 
As an example, consider the model in Sec. |j. With ay = (7r 2 /24)T 4 and choosing 
X = M x /T = f X i/(V2T) = a- 1 Xl ,weget 

g(x) = l + -^J ( X )-- 4 J (2x). (B.9) 

71* 71* 



Here Jo is from eq. QA.22 ). The numerical values following from the solution of eq. ( |B.2|) 



with this g(x) are shown in Fig. |^. Since fx = (a x /y / ay)/i = (4\/3//7r)(/i/T), it is seen 
from Fig. ^|that the asymptotic estimates in eqs. ( [B.7| ), ( |B.8| ) are reliable within ~ 10% 
at least up to /2 < 0.05, or fi/T< 0.02/. This corresponds to charges Q > 10 9 /" 4 . 
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